a=[0,1i];
b=a*(-0.5-0.88603i)+1i;
c=(a-1i)*(-0.5+0.88603i);
subplot(3,3,1);
plot([a,b,c]);
axis equal
for k=1:8
    subplot(3,3,k+1);
    a1=a/3;
    a=[a1,a1*(0.5+0.88603i)+1i/3,(a1-1i/3)*(0.5-0.88603i)+2i/3,a1+2i/3];
    b=a*(-0.5-0.88603i)+1i;
    c=(a-1i)*(-0.5+0.88603i);
    plot([a,b,c]);
    axis equal
end